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ABSTRACT 

Our previous theoretical study of the impact of an accreting envelope on the thermal state of an underlying white 
dwarf (WD) has yielded equilibrium core temperatures, classical nova ignition masses and thermal luminosities 
for WDs accreting at time averaged rates of (M) = 1 -1 1 - 1 0~ 8 M Q yr -1 . These (M) 's are appropriate to WDs in 
cataclysmic variables (CVs) of P or b < 7 hr, many of which accrete sporadically as Dwarf Novae. Approximately 
thirty nonmagnetic Dwarf Novae have been observed in quiescence, when the accretion rate is low enough for 
spectral detection of the WD photosphere, and a measurement of r e ff. We use our theoretical work to translate 
the measured T^'s into local time-averaged accretion rates, confirming the factor of ten drop in (M) predicted 
for CVs as they transit the period gap. For DN below the period gap, we show that if (M) is that given by 
gravitational radiation losses alone, then the WD masses are > 0.8M Q . An alternative conclusion is that the 
masses are closer to 0.6M Q and (M) is 3-4 times larger than that expected from gravitational radiation losses. 
In either case, it is very plausible that a subset of CVs with P or b < 2 hours will have Tiff's low enough for them 
to become non-radial pulsators, as discovered by van Zyl and collaborators in GW Lib. 
Subject headings: binaries: close — novae, cataclysmic variables- stars: dwarf novae — white dwarfs 



1. INTRODUCTION 

Cataclysmic variables (CVs; IWarnerl 119951) are formed 
when the low-mass stellar companion of a WD, exposed 
during a common envelope event, finally (on timescales of 
0.1-10 Gyr) fills the Roche lobe as a result of long term 
angular momentum losses. The WD will cool during this 
time ; a 0.20Afm He WD would have T c = 3.3 x 10 6 K at 4 
Gyr dAlthaus & Benvenutoll997h . where as a 0.6M^ C/O WD 
would have T c = 2.5 x 10 6 K in 4 Gyr (ISalaris et alJ 2000). 
These WDs have effective temperatures T e s « 4500 - 5000 
K just before mass transfer starts. However, once the Roche 
lobe is filled, the WD accretes material at (M) rs 10~ n - 
10" 8 M Q yr" 1 (e. g. Howell, Nelson, & Rappaport 2001) and 
can be reheated (Sio nll99ll) to higher r e ff's. 

Dwarf Novae (DN) are the subset of CVs with low time- 
averaged rates (M) < 10~ 9 M© yr" 1 and thermally unstable ac- 
cretion disks. The transfer of matter onto the WD occurs in 
outbursts that last a few days to a week once every month to 
year (or even longer in some systems). Most DN have orbital 
periods P 01 -b < 2 h ours, below th e "period gap", with fewer 
above the gap (see Shafterl [T992l) . During accretion disk qui- 
escence, the M onto the WD is often low enough that the sys- 
tem's UV (and sometimes optical) emission is dominated by 
light from the WD surface, allowing for a measurement of the 
WD r eff , nearly all of which exceed 10,000 K dSiorJI 19991) . 
Thus, the WD is hotter than expected for its age, providing 
evid ence of the thermal impact of prolo nged accretion o n the 
WD lSionll995l) . lTownsle v & Bildsten (2003) (hereafter^) 
have calculated T e s and its dependence on (M), the WD mass, 
M, and core temperature, T c . In this paper we now use that 
work to constrain these parameters from the measured r e ff 's. 

The gravitational energy released when a particle falls from 
a large distance to the stellar surface (GM/R) is deposited 



near the photosphere and is rapidly radiated away. This en- 
ergy does not penetrate inwards with the inflowing material, 
as the time it takes the fluid to move inward is much longer 
than the time it takes for heat to escape. This eliminates the 
outer boundary condition and instead points to the importance 
of energy release deep in the accreting H/He envelope due to 
both gravitational energy release and a low level of nuclear 
"simmering" (iTBl) . We begin in ^2] by reviewing our work 
and showing that the best constrained quantity from a mea- 
sured is (m) = {M)/A-kR 2 . At very low (M)'s, T efl also 
depends on T c . However, in this regime, we can calculate T c 
self consistently (TB), removing it from consideration. 

In |3]the (m)'s implied by the available measurements are 
presented and compared to CV evolutionary scenarios. A de- 
tailed comparison is presented for DN with P 01 -b < 2 hours, 
showing that the observed T^'s imply that either M > 0.6M Q 
or (M) is larger than implied by gravitational radiation losses 
alone. We close in 5|4]with a discussion of future work, espe- 
cially the seismology of accreting WDs. 

2. THE RELATION OF r EFF TO (M) 

The high quality UV spectra of quiescent DN from the STIS 
instrument on Hubb l e Spac e Telescope, e.g. iHowell et alJ 
(2002), Szkod vrtafl d2002c) and many of the references in 
the table appearing in Winter & Sion (2003), yield accurate 
measurements of T e s. However, the lack of accurate distance 
information prohibits the measurement of the WD radius, and 
hence mass, motivating the identification of a physical param- 
eter which is best constrained by T e s alone. The intent of this 
section is to demonstrate that, without knowledge of R or M, 
the best constrained parameter is the accretion rate per unit 
area, (m). 

Our previous work (TB) presented a detailed discussion of 
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FIG. 1 . — Left Panel: Predicted range of r cff (between lines) for 0.05M ien < M acc < 0.95M ign at M = O.6M (solid lines), l.OM (dashed lines), and 1.2M Q 
(dotted lines), Without knowledge of M, (m) is still fairly well constrained from T e ff. Right Panel: Dependence of T e ff range (as in left panel) on (M) for T c = T Cj eq 
(solid lines), 0.75r c . C q (dashed lines), and 0.5r c e q (dotted lines), at M = 0.6Mq . This observable does not depend strongly on T c except at the lowest (M) 's, and 
therefore is relatively insensitive to the CVs evolution. 



the impact of accretion on the thermal structure of a WD. For 
this paper we are primarily interested in the predictions for the 
surface luminosity, L{M, (M),T C ), for which we showed that 
the strongest dependence is L ~ (M)T C / /j,m p = AitR 2 <jsbT* s 
and thus oc (M)/4irR 2 = (m), even for different masses. 
The prime variance in this simple mapping is from the change 
in r e ff as the mass of the accumulated layer, M acc , increases 
between classical novae (CN). In discussing our theoretical 
results we use the T e s range for 0.05M; gn < M acc < 0.95M; gn , 
which represents where an observed WD is most likely to be 
found. We set Xj He = 0.001 throughout as the difference be- 
tween this and similar predictions for Xj Ue = 0.005 is less than 
the uncertainty due to the unknown M acc . Figure [2 shows 
the range of T e s traversed as a function of (m) for M = 0.6- 
1.2M©, showing that T & s provides a reasonable constraint on 
(m) even when M is not known. A T e ff = 15,000 K implies 
that 0.4 x 10" 3 < (m) < 1.3 x 10" 3 g cm" 2 s" 1 , whereas a 
r eff = 30,000 K implies that 1.0 x 10" 2 < (m) < 1.7 x 10" 2 

-2-1 

g cm s . 

The insensitivity of L to M is the strongest qualitative dif- 
ference between our calculation and Sion's (1995) estimate 
of L c ± 0.2( M)GM jR, made from the steady state models of 
llbenl (119821) . While Iben's (1982) models are steady state 
in the sense of L core = 0, they were for high accretion rates, 
(M) = 1.5 x 10~ 8 M© yr" 1 for the best discussed models, with 
one at (M) = 1.5 x 10~ 9 M Q yr" 1 , and use steadily burning 
shells. Such a steady state is inappropriate at the low (M)'s 
discussed here, where the burning is always unstable. Iben's 
(1982) models are also qualitatively very different from ours: 
in our models compressional energy released in the accreted 
H/He shell provides the outgoing luminosity and helps to es- 
tablish the equilibrium configuration of the core, whereas in 
Iben's (1982) models, the compressional heating term is en- 
tirely in the c ore, a contribution which is likely small (see 
Appendix A of TB). The results of TB are complementary to 
iGodon & Sionl (I2002T) . which focuses on the response of the 



WD to the short-timescale M variations during DN outbursts. 

The last parameter dependence to explore is T c , the WD 
core temperature. As discussed in lTBl DN below the period 
gap have adequate time to reach an equilibrium core temper- 
ature, r c eq , that depends on (M) and M. However, DN above 
the period gap have not been accreting long enough for T c to 
reach T C£q . Figure shows how the traversed T e g range de- 
pends on T c for a M = 0.6M Q WD. The curves show T e ff for 
T c = r Cj eq, 0.75r c eq and 0.5r c eq . Due to a strong core/envelope 
decoupling for (M) > 10" 10 M Q yr" 1 dH), the r eff range is 
nearly independent of T c . For lower (M)'s, the WD core tem- 
perature should be close to the equilibrium value, allowing us 
to use T c eq as representative when finding (m). 

3. INFERRING ACCRETION RATES FROM T EPP MEASUREMENTS 

Figure El s hows the (m)'s inferre d from the measured r e ff's 
tabulated in IWinter & Sionl (|2£)03), with the following mod- 
ifications: for WX Cet and VY Aqr we use the single- 
temperature fitted values, for CU Vel we correct the mis- 
quoted value, for AL Com we use a m easurement longer af- 
ter superou tburst (Szkodv et al. 2002d), and we add GW Lib 
JSzkodvetal.ll2002bl) and DW UMa ( ISzkodv et al.ll2002dl) . 
This observational (m)-P or b relation shows clear evidence for 
a drop in (m) below the period gap. 

The relationship between P 01 b and (At) is still a very ac- 
tive area of theoretical inquiry and one that we hope our work 
can ill uminate. The expec tations from "standard" CV evo- 
lution (iHowell et all 1200 It) for M = O.6M (solid line) and 
I.LMq (dashed line) are shown in Figure 2. In this disrupted 
magnetic braking scenario, (M) is set by magnetic braking 
above the period gap and by gravitational radiation below 
the period gap. Our deduced (wi)'s are lower than the ex- 
pected values above the period gap. It is important that we 
are inferring the long-term (M), averaged over the thermal 
time of the radiative (nondegenerate) layer ~ cpT c M n d/L « 
1O 4 ((M)/1O" 1O M yr" 1 )" 075 years for M = 0.8M Q (EH), so 
that this discrepancy cannot be due to a temporarily low M. 
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FIG. 2. — Values for the time averaged accretion rate per WD surface area, 
( m) = (M)/4nR 2 , derive d from the T e ff measurements in the table appearing 
in Winter & Sion 1 2003 ). The ranges indicated for each measurement are 
those allowed for 0.05M ig „ < M acc < 0.95M ign and 0.6M© < M < 1.2M© . 
The curves show the (m) predicted by Howell. Nelson. & RaDDarjort |20()1) 
forM = 0.6M© (solid line) and M = 1 . 1M© (dashed line). Patterson's (1984) 
deduced relation from CV observations is shown by the dotted line for M = 
0.6M© the dot-dashed line for 1 .0M© . The right hand scale gives (M)o,6, the 
corresponding accretion rate if R is that for M = 0.6M©. At the same (m), 

(M>i.o = 0.4(M)o. 6 . 

Although selection effects favor low (M) 's above the period 
gap, since such systems are more likely to have clean WD 
spectra, Patterson's (1984) estimates of (M) for the systems 
above the gap in Figure |2] are not systematically below his 
estimates for other CVs. The most recently improved cali- 
bration of the magnetic brakin g law, using spin-dow n of open 
cluster stars (Andronov. Pinsonnea ult. & Sillsll2003l) . yielded 
(M)'s above the period gap at least a factor of ten lower than 
Howell et al. (2001), falling below our inferences, so that 
braking in CVs must be enhanced over that responsible for 
the spin down of noninteracting low mass stars. 

In Figure 2, we also show Patterson's (1984) deduction 
from observations, (M) w 5.1 X 10" 10 (f or b/4 hr)M yr" 1 , for 
M = O.6M (dotted line) and l.OM (dot- dashed line). Our 
points are consistent with iPattersonl <fl9 84 ) above the period 
gap, within the uncertainty in his estimates. Below the gap, 
however, our measuremen ts are roughly a factor of 3 above 
his. The Patterson (1984) estimates also suffer from the ab- 
sence of reliable distances, but in a more direct way than ours, 
making systematic errors difficult to quantify. 

T he best quality data are f or CVs below the pe riod gap, 
and iHowell et all fc002l) and ISzkodv et all (l2002dl provide 
examples of these measurements with a discussion of how 
the uncertainty in the surface gravity, g, affects the re- 
sults. We display our predictions for the r e ff ranges along 
with observed values for P or \, < 2 hrs in Figure [5] The (M)- 
forb relation expecte d from gravitational ra diation losses for 
M = 0.6M Q is from Kolb & Baraffe ( 1999), and we use the 
same mass-radius relation for the donor to find (M)-P or \, for 
M= 1.0M Q . This gives (M) = 3.6 x 10" 11 and5.1 x 1O- U M 
yr" 1 for M = 0.6 and 1.0M Q respectively at P 01 -b = 1.5 hours. 
Due to a difference in the donor mass-radius r elation , these 
differ slightly from the values of Howell et al.l (120011) shown 
in FigureEl where (M) = 4.6 x 10" 11 and 6 x 10" n M Q yr -1 
for M = 0.6 and 1 . 1M Q at P orb = 1.5 hours. 

The measurement s show n by circles are again from the 
table in IWinter & Sionl |2003) (with the modifications dis- 
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FIG. 3 . — Comparison of our predicted ranges for T e ff with observed values 
for systems with P ol t < 2 hours. The filled areas indicate for each P OT t, the 
range of T e ff that a quiescent CV primary is expected to traverse between 
thermonuclear outbursts, using the (M) expected from angular momentum 
loss due to gravitational radiation iKolb & Baraffe 1999). Measurements are 
from the table in Winter & Sion 1 2003) except as noted in the text. Error bars 
indicate systematic errors due to the unknown WD mass, and represent the 
range of T e ff obtained with M = 0.3-0.9M© . All open points are subject to this 
same uncertainty. The filled points are systems which have well measured 
WD masses (Patterson 2001), in order of increasing P or ), and in units of M© 
the masses are 0.9 ±0.15, 0.82 ± 0.05, 0.61 ± 0.04, and 0.84 ± 0.09 fo r WZ 
Sge, OY Car, HT Cas, and Z Cha. The dia monds are magnetic CVs ISiorl 
rra^lGiinsicke et all200itlBelle et al f%M) . 

cussed earlier). The error bars indicate systematic errors due 
to the unknown WD mass, as the spectral fits cannot inde- 
pendently constrain T e ff and g, and represent the range of 
Jeff obtained for logg = 8 ± 0.5 (M = 0.3-0.9M©). All spec- 
tral measurements of these CV WD r e ff's are subject to this 
same uncertainty. The data with displayed error bars are the 
best measurements, and others have similar or larger u ncer- 
tainies. The diamonds show T e ff for magneti c systems JSionl 
U999tlGansicke et alJl200lUBelle et aT1l2003l) . The estimated 
masses for the two magnetic systems at 1 .9 hours are 0.5 and 
0.6M (n JSchwope et alJll993lTSchmidt, Stockman. & Grandil 
1983), for the lower and higher respectively, placing their 
measured T e ff very close to that expected from our work. 

When the WD mass, and thus radius, is known, (M) can 
be directly constrained. The filled points are systems which 
have relatively secure WD masses from eclipse timing (Pat- 
terson 2001), in order of increasing P or b and in units of M Q 
the masses are 0.9 ±0.15, 0.82 ± 0.05, 0.61 ±0.04, and 
0.84 ±0.09 for WZ Sge, OY Car, HT Cas, and Z Cha. The 
best r e ff measurement is that for WZ Sge JSion et alJI 19951) 
giving (M) = 6.4^} 4 x 1O" U M yr" 1 where the first er- 
rors represent the unconstrained value of M acc and the second 
the uncertain mass. If the gravitational radiation prediction 
of (M) l lKolb & Baraffdfl99^) is taken as a lower limit, then 
comparison of these CVs to the predicted r e ff's of Figure [5] 
yield a maximum WD mass. For example, none of the CVs 
denoted with filled circles can have M in excess of sa M , 
since a more massive WD would yield a higher T e ff than ob- 
served. 

If the current (M)-P or b relation predicted from gravitational 
radiation is correct, this comparison favors WD masses near 
0.9- 1 .0M Q . However, the systems with measured masses in- 
dicate another possible interpretation, that M « 0.85M Q for 
m any systems but ( M) is slightly greater than that predicted 
bv lKolb&Baraffel ( fl99^) . The fact that the lower bound of 
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measured values lies roughly parallel to what is expected for 
post-turnaround objects provides the exciting possibility that 
the turnaround is at roughly P 01 -b = 1.3 hours and the objects 
with the lowest T e g values, HV Vir and EG Cnc, are post- 
turnaround CVs. Under this interpretation, our work would 
indicate that the (M)'s for these post-turnaround objects are 
much higher than that expected from current modelling of 
evolution under the influence of gravitational radiation. Such 
an "extra " angu lar momentum loss has been mentioned by 
iPattersor] d200 11) as a way to understand the location of the 
period minimum. 

4. CONCLUSIONS 

We have used our theo retical work on accreting WDs 
jTownslev & Bildste nT2003l) to translate measurements of CV 
WD r e ff's into measurements of (m) = (M)/4ttR 2 , the accre- 
tion rate per unit WD surface area averaged over the thermal 
time of the radiative envelope ( > 1000 years). Since the pre- 
dictions of TB for L(M, (M) , T c ) are insensitive to T c (except 
where its value can be found self-consistently) we allowed 
for a range of WD masses, 0.6-1.2M Q , and accumulated en- 
velope masses, 0.05Mi en < M acc < 0.95Mi gn , to obtain a well- 
quantified uncertainty on (m) from the observed T e ff. We find 
evidence that above the period gap (P or b > 3 hours) DN ac- 
cret ion rates are slightly overestimated by CV evolution mod- 
els (iHowell et al.ll2001l) . but that angular momentum loss is 
quite enha nced compared to spin down of is olated low-mass 
stars ( Andronov. Pinsonneault. & Sills 2003). Below the pe- 
riod gap, P 01 b < 2 hours, we find that (M) is larger than that 
predicted by current models of gravitational radiation losses 
ilKolb & Baraffe- 119991) when M = 0.6M Q , indicating either 
larger M or higher (M) . 

It is well known that an isolated WD will pulsate 
when its T^ff is in t he approximate range 11000-12000 K 
( Berg eron et alJ fl995l) . While a difference in the outer at- 
mospheric composition, H/He in an accreting WD versus 
pure hydrogen in the isolated case, will shift this range, it 
is likely that a similar pulsation mechanism will be active 
in accreting WDs. Our calculations indicate that accreting 
WDs with M = O.6-1.OM should be near this range when 
(M) = fewxlO" n M Q yr" 1 . This (M) is typical of that ex- 
pected in CVs when accretion i s driven by gravitation al ra- 
diation losses, P 01 b < 2 hours jKolb&Baraffelll999f) . In 



fact one system, GW Lib, has been found which do es ex- 
hibit precisely this type of variability (van Zvl et aU fcOOO: 
Szkod v"et al.l 1200 2b!). Using the interior models developed 
in lTBl we are now undertaking a seismological study of these 
systems. This offers the tantalizing possibility of measuring 
the WD mass, spin and mass of the accumulated H/He layer. 

Since the minimum light of (M) < 10" 10 M Q yr" 1 
DN in quiescence is determined by the hot WD, we 
can now calculate just how deep the large-scale opti- 
cal surveys are probing into the predicted large popula- 
tion of CVs with very low mass co mpanions (< 0.1M Q ) 
jHowell. Rappapo~ & Politano| 119971) . A typical survey 
complete to V = 20 would find all 0.9M Q (0.6M Q ) WDs 
with T ef f > 8500 K (7000 K) that are within 200 pc 
(Bergeron, Wesemael, & Beauchamp 1995), and thus all DN 
with (AT) > 0.8 (1.0) x 10" n M Q yr -1 , including most DN be- 
low the period gap that have not yet reached the period mini- 
mum. Selection of the faint CVs is still a challenge; color- 
selected surveys utilize their unusual combinati on of a hot 
WD p lus a main sequence M star (see Townslev & Bildsten 
120021 for an earlier discussion and application to CVs in 
Galactic Globular Clusters). The second round of discover- 
ies of faint CVs from the Sloan Digital Sky Survey (SDSS), 
sensitive to ~20th magnitude, have just been announced 
dSzkodv et al.ll20 03). Based on t he initial set of 19 CV 's and 
their orbital period distribution, Szkod v"et alJ ( 1200 2 al) claim 
that a large fraction of the 400 eventual CV discoveries will 
be of t he low (M) variety below the period gap. Mar sh et all 
(2002) have also reported their discovery of three such sys- 
tems in the 2dF survey, sensitive to ^2 1st magnitude, and 
expect to have about 20 low (M) systems when the survey is 
complete. 
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